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Abstract. 

Several algorithms for tracking and for primary and secondary vertex reconstruction have been developed by the ATLAS 
collaboration following different approaches. This has allowed a thorough cross-check of the performances of the algorithms 
and of the reconstruction software. The results of the most recent studies on this topic are discussed and compared. 
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INTRODUCTION 

ATLAS is one of the four experiments at the Large Hadron Collider (LHC) that will produce pp collisions at a 
centre-of-mass energy of 14 TeV. There will be a bunch crossing every 25 ns, which imposes fast response from the 
detectors as well as the need to store data from each collision in on-detector pipelines until a decision of the first level 
of the trigger is taken. Furthermore there will be a very high density of charged tracks: at the nominal luminosity 
of 10 34 cm _2 sec _1 for each bunch crossing there will be about 200 charged tracks and about 15 vertex candidates. 
The reconstruction of tracks and of primary and secondary vertices at LHC will be a challenging task. The Inner 
Detector (ID) is the tracking detector of ATLAS and consists of three sub-detectors, whose performances as simulated 
by the software are extensively described in ([II 01- The outermost detector is the Transition Radiation Tracker (TRT) 
consisting of several layers of 4 mm straws in the barrel region (arranged in 3 layers of modules) and 14 Transition 
Radiation Tracker wheels in the endcap, providing about 30 hits per track and a resolution in the R(j> plane of 170 jXm. 
The Semiconductor Tracker (SCT) is located in the region between 25 cm to 50 cm in radius and consists of 4 layers 
of stereo silicon strips detectors in the barrel and 9 disks per side in the endcaps, achieving a resolution of about 17 
flm and 580 fltn in the R<p and Rz plane, respectively. The innermost is the pixel detector [3], consisting of 3 layers of 
silicon pixel detectors (at radius of 5.05, 8.85 and 12.25 cm) in the barrel and 3 disks in each of the endcaps, providing 
a track resolution of about 1 1 fim in R<p and 60 jim in Rz- To reconstruct the tracks, it is essential to take into account 
the effect of multiple scattering and of the energy loss in the material, which requires a precise knowledge of material 
in the detector. More complicated tracking algorithms are needed to consider the track resolution degradation at the 
edges of the ID due to the non-uniformity of the 2T magnetic field provided by the ATLAS solenoid in the ID region. 

TRACKING AND VERTEXING AT TRIGGER LEVEL 

The ATLAS trigger is subdivided into three different trigger selection layers that reduce the 1 GHz interaction rate to 
200 Hz. 

• The first level trigger (LVL1) is a hardware trigger with a 2.5 ms latency that brings the rate to 75kHz. It uses 
reduced granularity data of the calorimeter and the muon detectors and identifies geometrical regions of interest 
(Rol) in the detector. 

• The second level of the trigger (LVL2) uses the ID information since it processes in parallel the full information 
of all sub-detectors in the Rols defined by LVL1 . It has a latency of 10 ms bringing the rate below 2 kHz. Tracking 
and vertexing can be performed at LVL2. 

• Finally the Event Filter (EF) can use algorithms similar to the offline software having a 2s latency time. 

The performance of the tracking and vertexing algorithms implemented at ATLAS at LVL2 is described below. 



Tracking at LVL2 



The tracking algorithm forms track seeds by fitting with a straight line pairs of space points in the pixel innermost 
layer (B layer) and in the second logical layer (in a given Rol). The tracks are extrapolated back to the beam line and an 
Impact Parameter (IP) is obtained for each of them. The track is retained if the IP in the transverse plane with respect to 
the beam axis is small. The z coordinate of the primary vertex is then obtained as the maximum of the histogram filled 
with the z intersection of the seeds with the beam line. A third space point is extracted in modules situated in positions 
where the hits may lay based on the track extrapolation. After having removed the ambiguities due to overlapping 
space points in triplets using the extrapolation quality, the triplets are fitted and identified with tracks. The efficiency 
for tracks in jets is 80-90% depending on the luminosity and the event topology and 95% for single electrons. 



^-tagging at LVL2 

The b-jet selection is performed by using the transverse impact parameter do significance, defined as S = do/s(do), 
where s(do) is the error on do and its dependence on pj is obtained from the simulation. 

A secondary vertex algorithm similar to offline but faster is used. The Z?-jet estimator uses a likelihood ratio given 
by the product of the ratios of the probability densities for each track to come from a b-jet or a light-jet, 

w = n/ b(Si) 



MSiY 



where i runs over all reconstructed tracks in the jet. 

The final discriminative variable is defined as X = W/(l+W). The processing time for fe-jets is less than 2 ms. The 
rejection for light quarks is 25% and 15% for a Z?-jet efficiency of 50% and 60% respectively. 



OFFLINE TRACKING ALGORITHMS 



Several offline tracking algorithms have been developed at ATLAS: xKalman and iPatRec and only very recently the 
NewTracking algorithm. 

• The xKalman first searches for tracks in the TRT using fast histogramming of straw hits, then it extrapolates back 
to SCT and pixels and fits the tracks using a Kalman filter that associates clusters to tracks assuming that the 
noise and all material effects and measurements are Gaussian. Finally, the improved tracks are extrapolated back 
into the TRT in a narrow region around the extrapolated trajectory, retaining all hits in that region. 

• the iPatRec algorithm forms track-candidates in SCT and pixel detectors using space-point combinatorials subject 
to criteria on maximum curvature and crude vertex region projectivity. A global % 2 fitter is used to fit tracks and 
associate clusters. Only good tracks are retained for extrapolation in TRT, where TRT hits are added. To limit the 
contamination from high occupancy, tight cuts are applied on the straw residuals. 

• The NewTracking algorithm is mainly a reorganization of the tracking code. At present it is largely based on 
xKalman, but in the future it will use also tools from iPatRec with the aim of obtaining optimised performances. 
Furthermore it uses a better detector geometry description. 

Table Q] shows the results of a performance comparison done using tt events. As can be seen xKalman and iPatRec 
have comparable performances, slightly better than the newly developed NewTracking, whose performances are going 
to improve and are already better than the others at the time of writing. 

Figures Q] a) and b) show the resolution on the impact parameter in the transverse and longitudinal planes, respec- 
tively. A sample of WH events with m//= 400 GeV /c 2 has been used for the study, where the W decays to ji and H 
to mm. Figure [2] shows the momentum resolution versus pj for single muons averaged over all |tj | and for rj=0. 



PRIMARY VERTEX RECONSTRUCTION 



Given the large multiplicity of tracks at each bunch crossing (several hundreds as discussed before) the vertex 
reconstruction must be fast and robust. The input consists of the 3 dimentional trajectories and error matrices of 



TABLE 1. Performance of the ATLAS offline tracking algorithms ob- 
tained using tt events. The result of the newly developed NewTracking are 
very preliminary and already superseeded at the time of writing. 



xKalman iPatRec NewTracking 



Multiplicity (p>l GeV) 16.69 17.06 16.88 

Barrel Track eff/fake rate 99%/0.6% 99%/0.7% 96%/2.5% 

Transition eff/fake rate 98%/0.6% 98%/0.5% 96%/3.6% 

Forward eff/fake rate 98%/0.3% 99%/1.3% 95%/2.7% 
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FIGURE 1. Figure a) and b) show the transverse and longitudinal impact parameter resolution for a sample of WH events with 
niH= 400 GeV/c , where the W is decaying to /ivy, and the H to uu. 

the tracks. Quality requirements on tracks are applied: typically pr > 1 GeV/c, Idol < 0.25 mm, Izol < 150 mm and % 2 
per track <0.5. 

The approximate primary vertex position in z is found using a sliding window of 0.7 cm that is moved along the 
whole interaction region. The window with largest number of tracks, weighted with pj is chosen. The <z> position of 
the vertex is given by the mean of all the tracks in that window. Tracks belonging to the primary vertex are taken away 
and the procedure is iterated to get other (pile-up) vertices. 

All tracks at ± 5 mm in z and ± 1 mm in the transverse plane are accepted as coming from primary vertex. At this 
point the vertex fitting is performed using a Billoir method. In the fitting procedure the outliers are removed, that is if 
the % 2 obtained when adding a track is too high, the track is rejected and the fit is recalculated, 

There are two different implementations of this method which are basically using the same strategy: VxPrimary and 
VKalVrt. 

A third primary vertex fitter has been developed, the Adaptive Vertex Fitter, that solves the problem of outlier tracks 
that spoil the fit, not by discarding them, but by down-weighting them. It minimises instead than residuals, the sum of 
squared residuals weighted with the % 2 . Table|2]shows a comparison of the performances of the primary vertex finder 
algorithms, which give very similar results by following slightly different approaches. 



ATLAS TAGGING 



Basic ingredients for the ^-tagging are the long lifetime of the fo-hadron jets, the high multiplicity Z?-jets and the 
impact parameters of the tracks. Additionally the presence of a secondary vertex can increase the rejection power of 
the b selection. The /?-tagging algorithms developed by the ATLAS collaboration are based on such properties of the 
fe-jet [4]. In addition ^-tagging methods based on the tagging of a soft lepton (electron or muon) resulting from the 
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TABLE 2. Primary vertex resolution in x 
and z for the ATLAS primary vertex finders. A 
sample of simulated WH events with m#=l20 
GeV/c 2 and W — > flVn, H —t bb has been used 
for the study. 







z((lm) 


VxPrimary 


12.6 ±0.1 


50.0 ± 0.5 


AVF 


11.07 ±0.09 


46.76 ± 0.05 


VKalVrt 


11.07±0.09 


45.43 ± 0.05 



decay of a ^-flavoured hadron are being developed. 

The standard fo-tagging identification is based on the 2-dimentional impact parameter of a track in the transverse 
plane. The method adopted by ATLAS uses the normalised significance S — do/s(do) for each track and compares it 
to predefined calibration probability density functions for the b and light quark hypothesis to obtain the probabilities 
b(S) and u(S). The following discriminating variable (IP2D) for each jet is built by summing over all the tracks in a 

"' ft K(S,) 

The discriminative variable consists of logarithms of the density functions ratio, allowing to combine two variables 
easily. The ^-tagging performance can be improved by adding the longitudinal and the transverse significance (IP3D): 

w = Pb{Sd ,Szo) 

Pu (Sd J Sj ) 

Figures[3]a) and b) show the distribution of the transverse and longitudinal impact parameter significance for b and 
light quarks, respectively. 

In addition the presence of a secondary vertex is searched for. Only tracks passing the following quality cuts are 
used: pj > 1 GeV/c, |tj| < 2.5, |iio| < 1 mm > kol < 1-5 mm, number of pixel hits in the B layer greater than and 




number of pixel hits larger than 1, while more than 6 SCT hits should be detected. The selected tracks are used to 
search for good 2 track vertices in the jet. At this point vertices that are due to decays of Kaons or Lamdba baryons or 
to interactions with the beam pipe material, the pixel layers or to photon conversions (VO's) are removed. A common 
(inclusive) vertex is formed for the remaining tracks. 

Once the secondary vertex is reconstructed some additional discriminating variables can be considered: 

1. the number of good two track vertices in the jet, 

2. invariant mass of all particles coming from the secondary vertex mass, 

3. the ratio of the energy of all particles coming from the secondary vertex and the total energy of the jet. 

The variables have been chosen to be independent from tracks impact parameters in order to have a real gain on b- 
tagging performance. The secondary vertex discriminating variable obtained from the probability density functions 
that parametrize the distributions of these variables can be combined with the track impact parameter based ^-tagging 
procedure. Figure|4]a), b and c) show the distribution of the secondary vertex variables for b and light jets. 

A jet is labeled as a fo-jet if there is a fo-quark within a cone of radius 0.3 around the jet axis, the efficiency for fo-jets 
(£/,) is defined as the ratio of the number of jets and the number of jets labeled as b with pj > 15 GeV/c and | tj | < 2.5. 
The light jet rejection is defined as R u = 1/ £„. 

Overlapping jets cause mislabeling and furthermore the jet isolation is very dependent on the type of physics process 
taken into account. A purification of the jets is done to factorize those effects from pure fo-tagging issues: light jets 
are not taken into account if there is a fo/c/quark/hadron within a cone of radius 0.8 around the jet axis. Table|3]shows 
the rejection for light jets for a Z?-tag efficiency of 50% and 60% for the IP2D and IP3D impact parameter and for 
the combined IP3D and SV1 discriminating variables. Figure [4] shows the combined SV1+IP3D weight. The results 
are obtained using a sample of WH-» jj.v^uu events with m// =120 GeV/c 2 . Tracks have been reconstructed with 
the xKalman algorithm and two different primary vertex finders have been compared: the adaptive vertex fitter and 
VkalVrt. In the last two rows of the table, results are shown when rejecting bad tracks recognised as coming from VO's 
or interactions with the beam pipe or detector material. In this study the jets were reconstructed around the primary 
quark directions with a cone size of AR= 0.4. 

An evaluation of the /?-tagging performance using 190K of tt events has shown similar results: the rejection for 
light quarks obtained using SV1+IP3D for a 50% efficiency for fr-quarks is 858 ±42.9 while for a 60% efficiency is 
259 ±7.8. 
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FIGURE 4. Figure a), b), c) show the distribution of the secondary vertex variables for b and light jets. Figure d) shows the 
weight obtained by combining SV1 and IP3D. 

TABLE 3. b-tagging rejection for an efficiency of 50% and 60% for the IP2D, IP3D and SV1 taggers, 
using the AVF and VkalVrt primary vertex fitters. A sample of WH^ pLV^uu events with m# = 120 
GeV/c 2 events has been used for the study. The last two rows show the results when rejecting tracks 
coming from VO's. 



IP2D IP3D IP3D+SV1 



efficiency 


50% 


60% 


50% 


60% 


50% 


60% 


Rej. VKalVrt 


135±9 


55 ±2 


214 ± 18 


75 ±4 


609 ± 86 


157 ± 11 


Rej. AVF 


130 ±9 


52 ±2 


205 ± 17 


73 ±4 


612±87± 


147 ±10 


Rej. no VO's + VkalVrt 


206±17 


69 ±3 


339 ±35 


101 ±6 


815 ± 134 


192 ± 15 


Rej. no VO's + AVF 


199 ± 16 


66 ± 3 


327 ± 34 


98 ±6 


794 ±129 


164 ± 12 



CONCLUSIONS 

There has been a lot of work and improvements on the tracking and vertexing algorithms in the ATLAS collaboration 
in the past few years. Different approaches have been followed in parallel to develop the tracking and vertexing 
algorithms, giving comparable results. Cosmic events have successfully been reconstructed with the SCT and TRT 
barrel on the surface in the spring of this year. 
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